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Abstract 

QCD sum rules for D and B mesons embedded in cold nuclear matter are evaluated. 
We quantify the mass splitting of D — D and B — B mesons as a function of the nuclear 
matter density; extrapolated to saturation density it is in the order of 60 and 130 MeV, 
respectively, driven essentially by the condensates (q^q), (q^gaWq) and (qq}. The genuine 
chiral condensate (qq), amplified by heavy-quark masses, enters the Borel transformed 
sum rules for the mass splitting beyond linear density dependence. Including strange 
quark condensates reveals a numerically smaller and opposite effect for the D s — D s mass 
splitting. 



I. INTRODUCTION 



QCD sum rules offer a link from hadronic properties, encoded in spectral func- 
tions, to QCD related quantities, like condensates, in the non-perturbative domain. 
A particularly valuable aspect of QCD sum rules is, therefore, the possibility to 
predict in-medium modifications of hadrons, supposed the density and temperature 
dependence of the relevant condensates is known. Taking the attitude that this 
is the case, one arrives at testable predictions for changes of hadronic properties 
in an ambient strongly interacting medium. There is a vast amount of literature 
on the in-medium changes of light vector mesons, cf. [IJEJEIIIIEJEIEIEIIH] an< ^ 
further references therein. Vector mesons are of interest as their spectral functions 
determine, e.g., the dilepton emissivity of hot and compressed nuclear matter. Via 
the direct decays V — > where V stands for a vector meson and for a 

dilepton, the spectral distribution of V can be probed experimentally. Accordingly, 
heavy-ion experiments are often accompanied by special devices for measurements 
of l + l~ = e + e~ or /i + yU~. Addressed questions concern in particular signals for chiral 
restoration [8]. Clearly, besides the QCD sum rules, also purely hadronic models 
have been employed to understand the behavior of vector mesons in nuclear matter, 
cf. [U [U |9] for examples. 

Such hadronic models are also used in the strangeness sector [10J. Here, the 
distinct behavior of kaons and anti-kaons attracted much attention, cf. [TT] for ex- 
perimental aspects. The upcoming accelerator complex FAIR at GSI/Darmstadt 
offers the opportunity to extend the experimental studies into the charm sector. 
The CBM collaboration [12] intends to study the near-threshold production of D 
and J /if) mesons in heavy- ion-collisions, while the PANDA collaboration [13] will fo- 
cus on charm spectroscopy, as well as on charmed mesons produced by anti-proton 
annihilation in nuclei. In the CBM experiments, charm degrees of freedom will serve 
as probes of nuclear matter at the maximum compression achievable in the labo- 
ratory, at moderate temperatures. Despite of this interest in D mesons and their 
behavior in nuclear matter, the literature on in-medium D mesons is fairly scarce. 
While there is a variety of calculations within a hadronic basis, e.g. [T3] IT5] [16], or 
within the quark-meson coupling model, e.g. [T7], the use of QCD sum rules is fairly 
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seldom [18j [T9l [20] . In contrast, the treatment of vacuum D (and D s ) ground states 
is performed in a concise manner (2U [221 [23] • 

The aim of the present paper is the re-evaluation of the QCD sum rules for 
D and D mesons in cold nuclear matter and an extension to B and B mesons 
as well. Even for the operator product expansion (OPE) up to mass dimension 
5, there are conflicting results in the literature concerning the open charm sector 
PH [HI ESB 1231 [21 123 125] • While in jH] only the even part of the in-medium OPE 
up to mass dimension 4 has been used, we present here the even as well as the odd 
in-medium OPE up to mass dimension 5. Moreover, a term oc (qga&q), i.e. the 
lowest-order quark-gluon condensate, can be found in the literature with various 
factors and signs already for the vacuum. As the subtle D — D mass splitting is of 
paramount experimental interest, a safe basis is mandatory. 

Our article is organized as follows. Section [IT] contains the QCD sum rules for- 



malism for D and D mesons. The spectral functions are discussed in section III 



The numerical evaluation for D, D and B, B mesons is executed in sections IV and 



[V] The discussion and summary can be found in section VI 



II. QCD SUM RULES 

The basic quantity to be evaluated is the two-point function 

n(g) = i I d 4 xe^(Q\T [j(x)jt(0)j |0) (2.1) 

as the Fourier transform of the expectation value of the time-ordered product of 
the currents j(x) and j^(0); the state has properties J^\Q) = Eq\H), (0(0) = 
1, a\Q) 0. M' is the full Hamiltonian of the theory, a an arbitrary annihilation 
operator, and the field operators are taken in the Heisenberg picture. Splitting up 
n(g ,<f) into an even (e) and an odd (o) part according to U(q ,q) = U e (qQ,q) + 
%n°(gg,«f) with 

n e (g , q) = \ (n(g , q) + n(-g , q)) = II e (-g , q) , (2.2a) 
n°(go, q) = ^~ (n( 9 o, q) - n(-g , q )) = n°(-<7o, q) , (2.2b) 
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one arrives at the TV-fold subtracted dispersion relations in the complex q plane 
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with AIl(s) = Imll(s). From Eq. (2.3) we see that H e '° depend on The Borel 
transformed sum rules are 
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where the subscript OPE denotes the operator product expansion of 
(0|T [j(a;)j%)] = T.o°o{x - y)(tt\0\Q) with QCD condensates (tt\0\Q) 
and Wilson coefficients Co- We are interested in the low- lying strength en- 
coded in J s " ds AUse~ s2 / M2 and f s ° ds AHe~ s2 ^ M2 , while the continuum parts 
(l-L + I4 ) ds Anse" s2 / M2 and ( + f s +°°) ds AUe~ s2 / M2 will be merged into 
the perturbative OPE part 14^+ (s) (see below) according to the semi-local duality 
hypothesis; Sq are the corresponding continuum thresholds; M is the Borel mass. 

Employing the current operator j D + = id^ 5 c (and j D - = j^+(x) = icq^d for the 
antiparticle), we obtain for the OPE side up to mass dimension 5, in the rest frame 
of nuclear matter v = (1,0) (v stands for the medium four- velocity) , in the limit 
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and sufficiently large charm-quark pole mass m c , 
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where a s = g 2 /47r. (Analog relations hold for ] D o(x) = iujsc with j^o(x) = j 



icj^u.) The calculational details are documented in [27]. While the perturbative 
spectral function Imll^+(s) (see [23l [24] for an explicit representation in terms of 
the pole mass) is known for a long time, discrepancies especially for Wilson coeffi- 
cients of medium specific condensates exist. An important intermediate step is the 
careful consideration of the operator mixing, which occurs due to the introduction 
of non-normal ordered condensates and the corresponding cancellation of infrared 
divergent terms cx m~ 2 and log m q (m g is the light-quark mass) at zero and non-zero 



densities • This is not to be confused with the operator mixing within renormal- 
ization group methods. In vacuum our expression differs from [TB] in the prefactor 
of ((a s /ii)G 2 ); [23] reports an opposite sign; [19] finds the same result. For the 
medium case [19] does not give explicit results, while terms oc (dd), oc ((a s /ii)G 2 ), 
oc ((a s /n)((vG) 2 /v 2 — G 2 /4)) have different prefactors compared to [18]. Higher 
order terms are partially considered in [19] and are found to be numerically not 
important. 

We stress the occurrence of the term m c {dd). In the pure light quark sector, say 
for vector mesons, it would read rrid{dd), i.e., the small down-quark mass strongly 
suppresses the numerical impact of the chiral condensate (dd). In fact, only within 
the doubtful factorization of four-quark condensates into the squared chiral conden- 
sate it would become important [7J. Here, the large charm-quark mass acts as an 
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amplifier of the genuine chiral condensate entering the QCD sum rules for the D A 
meson. 

III. PARAMETERIZING THE SPECTRAL FUNCTION 



Especially in vacuum the spectral strength of the iso-scalar-vector excitation ex- 
hibits a well-defined sharp peak (the u meson) and a well-separated flat continuum. 
Assuming the same features for the u meson in a medium gives rise to the often 
exploited "pole + continuum" ansatz. One way to avoid partially such a strong 
assumption is to introduce certain moments of the spectral function, thus replacing 
the assumed pole mass by a centroid of the distribution [7J 128] . 

For D mesons the sum rule includes an integral which arises from the dispersion 



relation over positive and negative energies, see Eq. (2.4). Similar to baryons 



one may try to suppress the antiparticle contribution corresponding here to D~ . 
This, however, is not completely possible [6]. Nevertheless, one can identify with 
the ansatz AIl(s) = nF + 5(s — m + ) — -kF- 5(s + m„), motivated by the Lehmann- 
representation of the correlation function, the meaning of the even and odd sum 



rules (2.4) with (2.5): 



ds s AIIe- s2 / M2 = m + F + e- m2 +/ M2 + m„F_e- m -/ M \ 
"° + ds AITe- s2 / M2 = F+e- m l' M2 - F.e"™ 2 ^ . 



(3.1a) 
(3.1b) 



With the decomposition m± = m ± Am and F± = F ± AF the leading order terms 
of an expansion in Am for the first and second lines become cx Fme~ m2 ^ M2 and 
oc (AF — 2AmF|j)e" m2 / M2 meaning that (3.1a) is related to the average D + D 



properties, while (3.1b) refers to the D — D splitting. If one assumes for the moment 



being m± and F± to be independent of the Borel mass M, (3.1 ) can be rewritten as 
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where a prime denotes the derivative w.r.t. 1/M 2 . In order to gain further insight 
into the dependencies of Am and m on the different OPE contributions, we expand 



(3.2) up to first order in the density n employing e(n) pa e(0) + nde/dn\ n=0 and 
o(n) pa ndo/dn\ n= o, since o(0) must vanish to reproduce the vacuum sum rules 
where Am(n = 0) = holds. We remark that these expansions are exact for a 
linear density dependence of the condensates and if Sq = ((sq ) 2 + (sq ) 2 )/2 as well 
as Asq = ((sq ) 2 — (sq) 2 )/2 are density independent and the Borel mass M is kept 
fixed. This implies As 2 , = for all densities, because otherwise o(0) = cannot be 
fulfilled. For small densities we get accordingly 

1 #Le'(0) -e(0)^L 
Am(n) Pa ^ e (0V» n > ( 3 - 3a ) 



, , . : e'(0) , 1 / e(0) £| e'(0)-e(0)g| 
m{n) pa \ y-- + -W — — y — — ^n, (3.3b) 



which can be written as 



e(0) 2 V e'(0) e(0) 2 



l^lo^(°) + f lo. 



Amfn) « — _onjo LjI (3. 4a ) 

V ; 2 e(0) ' V ; 
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V ; V ; 2m(0) e(0) V ; 

Eq. ( |3.2 ) and the approximations in (3.4) offer a transparent interpretation. In 
vacuum (n = 0), there is no mass splitting, of course; the mass parameter m(0) is 
determined by the even part of the OPE. In first order of n, the mass splitting Am 
depends on both the even and odd parts of the OPE, whereas only the even part 
of the OPE determines the mass parameter m, having the meaning of the centroid 
of the doublet D + , D~ . If one is only interested in the mass shift of the iso-doublet 
as a whole, for small densities it is sufficient to consider the even OPE part alone, 
as was done in [18j. However, for the mass splitting the odd part of the OPE is of 
paramount importance. In particular, it is the density dependence of the odd part of 
the OPE alone which drives the mass splitting in first order of n. Interestingly, the 
density dependent part of the chiral condensate, which belongs to the even part of 
the OPE, enters the mass splitting in order n 2 . The chiral condensate comes about 
in the combination m c (dd). The large charm mass amplifies the numerical impact, 
as stressed above. 
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Furthermore, up to order n, only sl(n = 0) and 
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seen from the derivatives needed to calculate m and Am from (3.3): 
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While (3.4) suggests that one can independently adjust m(0) to the respective 



vacuum value, Eq. (3.5) evidences that further vacuum parameters (such as M, 

) enter the density dependence and have to be chosen 
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consistently to the vacuum mass. That means, one has to evaluate the complete 
sum rule, including consistently the vacuum limit. 



We remark that (3.2) or (3.4) are a consequence of using a pole-ansatz for the 



first excitation. The OPE and the special form of the continuum contribution to 
the spectral integral are encoded in e and o. Likewise, the arguments following 
( 3.2[ ) merely use o(0) = 0. The last point must always be fulfilled in any sum rule 
and/or dispersion relation, because at zero density, the current-current correlation 



function (2.1 ) only depends on q and, hence, the odd part (2.2b) vanishes. This can 



also be confirmed directly from (3.1b), where Sq = Sq , due to particle anti-particle 
symmetry, and AIl(s) = AIl(s 2 ), meaning that the spectral density in vacuum 
merely depends on the squared energy, on account for o(0) = 0. 

To arrive at a more general result, one may seek for a relation of m± to certain 
normalized moments of H(s) (or ratios thereof) independent of a special ansatz, as 
can be done in the case of vector mesons [71 [2H] • In this spirit one would be tempted 
to define f*° ds sAIie- 3 ' 2 ' M ' 2 -> m+F +e - m + /M2 and f*° dsAUe~ s2 / M2 -> F + e~ m + /M2 
and analogously for m_ and _F_ . However, such a separation of positive and negative 
frequency parts leads to multiple but different expressions for m± which can be 
fulfilled consistently only for special cases of II(s), as for the above pole ansatz. 
(This can be seen by combining these relations with derivatives according to M~ 2 .) 



Therefore, one is left with either the somewhat vague statement that (3.1) refers to 



S 



D + D and D — D properties or one has to employ another explicit ansatz for the 
function U(s). 



Alternatively, one can define moments which correspond to the integrals in (3.1) 
S n (M) = f ° dss n ATl(s)e- s2/M2 . (3.6) 

Js~ 

The odd and even OPE, o = So(M) and e = — S\(M), and their derivatives with 



respect to M 2 , o' = —Ss(M) and e' = S^M), can then be related via (3.2) to 



these moments. Thereby, new quantities Am and m may be defined which encode 
the combined mass-width properties of the particles under consideration: 



1 S1S2 — SqS[ 
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and m 2 = Am + m + m_. For the above pole ansatz, these quantities become 



Am = Am and m = m, i.e., they allow for an interpretation as mass splitting and 
mass centroid. The relations (3.6) and (3.7) avoid the use of a special ansatz of the 
spectral function, but prevent a direct physical and obvious interpretation. 

IV. EVALUATION FOR D AND D MESONS 

We proceed with the above pole ansatz and evaluate the behavior of m± having 
in mind that these parameters characterize the combined D, D spectral functions, 
but need not necessarily describe the pole positions in general. According to the 
above defined current operators, D stands either for D + or D° and D for D~ or D° . 



Because dm±/dM = has been used to derive (3.2) we have to look for the 
extrema of m±(M). Furthermore, in order to solve consistently the system of equa- 
tions defined by ( 3.1[ ), the values taken for m± must be fixed at the same Borel 



mass M. Therefore, we evaluate the sum rules using two threshold parameters 
(sq ) 2 = Sq ± Asq and demand that the minima of the respective Borel curves 
m + (M) and m_(M) must be at a common Borel mass M. Hence, the thresholds 



are prescribed and offer the possibility to give a consistent solution to (3.1). 

Analog to the analysis in [18], we chose the threshold parameter Sq = 6.0 GeV 2 , 
which approximately reproduces the vacuum case. At zero density we obtain for m± 



TABLE I: List of employed condensate parameters. A discussion of these numerical values 
can be found in }29\j : further remarks on (q^gaWq) are given in [19] ' . For the strong coupling 
we utilize a s = 4ir/ ((11 — 2Nf/3) ln(n 2 / Aq CD )) with fi being the renormalization scale, 
taken to be of the order of the largest quark mass in the system, and Nf being the number 
of quark flavors with mass smaller than /u; Aq CD = 0.25 GeV 2 is the dimensional QCD 
parameter. The employed quark pole masses are m c = 1.5 GeV and rrib = 4.7 GeV ]25^ . 
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(-0.245 GeV) 3 


45/lln 
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-0.65 GeVn 
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-0.05 GeVn 


(qUDoq) 





0.18 GeVn 


(q [Dl - \go<S] q) 





-0.3 GeV 2 n 


(d ] Dlq) 





-0.0035 GeV 2 n 







0.33 GeV 2 n 



a value of 1.863 GeV, representing a reasonable reproduction of the experimental 
value of the D mass. The employed condensate values are listed in Tab. [TJ 

The density dependence of the mass splitting parameter Am and the D + D 
doublet mass centroid m are exhibited in Fig. [T] as a function of the density. We 
observe an almost linear behavior of the mass splitting with increasing density. At 
n = 0.15 fm -3 a mass splitting of 2Am w —60 MeV is obtained. The mass splitting 
has negative values, i.e. m_ > m + or mg > rrirj in line with previous estimates 
in [12]. For the mass centroid m our result differs from the one in [TS], where a 
mass shift of the order of —50 MeV is obtained, while we find about +45 MeV. At 
n = 0.15 fm~ 3 the splitting of the threshold parameters is As 2 , w —0.3 GeV 2 for the 
used set of parameters, and the minima of the Borel curves are located at M ~ 0.95 
GeV and are slightly shifted upwards with increasing density. 

While the mass splitting is fairly robust, we find a sensitivity of the centroid mass 
shift under variation of the continuum threshold parameter S 2 ,. The above reported 
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0.0 0.05 0.10 0.15 ' 0.0 0.05 0.10 0.15 



n [fm -3 ] n [fm -3 ] 

FIG. 1: Mass splitting parameter Am (left) and mass centroid m (right) for D, D mesons 
for density independent threshold (solid line) and a density dependent threshold Sq(u) = 
■Sq(O) ± n/no GeV 2 , where the dotted (dashed) curve is for the positive (negative) sign. 
Note that the mass splitting in the iso-doublet is 2 Am. 

value of the mass centroid changes towards zero when lowering Sq. In Fig. [I] we 
therefore also use a density dependent prescription for the threshold s§(n) = Sq(0) ± 
n/riQ GeV 2 , where no = 0.15 fm -3 is the nuclear saturation density; ±l/no GeV 2 
corresponds to the first Taylor coefficient ds^/dn^). This simple choice enables us 
to identify the uncertainties which might emerge due to the introduction of a density 
independent threshold. As can be seen, the average mass shift may change in sign. 
In contrast, the result for Am shows only a weak dependence on s"q. 

At this point a comment concerning the sign of (q^ga&q) is order. If one would 
use (q^ga&q) = —0.33 GeV 2 n instead (this option is also discussed in [29], (q^D^q) 
would acquire a value of —0.0585 GeV 2 n accordingly) one would get a much larger 
mass splitting of about —180 MeV, which is far beyond the estimates obtained in 
[TU [El HE]. Hence, we favor the positive sign of (q^gaWq) as advocated in [19], 
too. Clearly, further correlators should be studied to investigate the role of the 
condensate (q^ga&q). 

We emphasize the special evaluation strategy employed so far. Other possibilities 
are, e.g., variation of Sq and As 2 so that m±(M) develop a section of maximum 
flatness. Interestingly, this method leads to a rather low threshold Sq ~ 4 GeV 2 and 
a low vacuum mass of about m ps 1.6 GeV. In contrast, averaging over the Borel 
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FIG. 2: Borel curves m±{M) for the D meson for two densities and two values of the 
chiral condensate. Left panel: n = 0, right panel: n = 0.15/m™ 3 . (Solid curves: chi- 
ral condensate from Tab. [7| dotted curves: doubling the chiral condensate (left panel) or 
doubling the density dependent part of the chiral condensate (right panel); lower (upper) 
curves in the right panel are for m + (m^), while m + = m_ for the vacuum case in the 
left panel.). 

curves in the interval [0.9Mo, 1.2Mq], around the minimum Mo, we find the values 
for the mass splitting Am m —40 MeV and the average mass shift to be of the same 
order as quoted above, whereas the absolute value of the vacuum mass becomes 
m = 1.877 GeV. 

Let us now further consider the impact of various condensates. The result for 
the mass splitting Am strongly depends on the quark density (q'q), whose density 
dependence is uniquely fixed. The odd mixed quark-gluon condensate (q^gaGq) and 
the chiral condensate (qq) are the next influential ones for the mass splitting. The 
density dependent part of the chiral condensate enters in order 0(n 2 ) gaining its 
influence from the heavy quark mass amplification factor. The influence of the chiral 
condensate is illustrated in Fig. [2j In a strictly linearized sum rule evaluation, the 
density dependent part of m c {qq) would be omitted for the mass splitting. However, 
numerically the influence of the chiral condensate is of the same order as (but still 
smaller than) the above discussed condensate (q^ga&q), which enters the odd part 
of the OPE. As expected, the density dependence of the mass centroid is basically 
determined by the even part of the OPE. The density dependent parts of the other 
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even condensates are of minor importance for the mass splitting. The shift of the 
centroid's mass is anyhow fragile. 

Within the given formulation and with the first evaluation strategy, one may 
also consider D s and D s mesons with the replacements m q — > m s , (qq) — > (ss) = 
0.8(m)vac + y(qq)med, (qg^q) -> (sgo&s) = 0.8 GeV 2 (ss), (q^q) -> (sU) = 0, 
(qUD q) -> (sUDqs) = 0.018 GeVn, (q [D% - |</a^] g) -> (s [D 2 - \g<&] s) = 
y{q[D*-lgo9]q), (qW 2 q) - (s^o s ) = 2/(^o9>> (?V^<7> - (^V^) = 
y(q* ga&q) . The anomalous strangeness content of the nucleon is varied as < y < 
0.5 [31J; lattice calculations, for example, point to y = 0.36 [32]. The results are 
exhibited in Fig. [3] At n = 0.15 fm -3 and y = 0.5 we observe a mass splitting of 
2Am pa +25 MeV and a shift of the mass centroid of about +30 MeV. The splitting 
of the thresholds becomes As 2 Q pa 0.83 GeV 2 , and the minima of the Borel curves 
are located at M ~ 0.89 GeV and slightly shifted upwards with increasing density. 
The main reason for the positive sign of the mass splitting is the vanishing strange 
quark net density (s^s). The mass splitting acquires positive values for (s^s) < 0.4 n 
(at y = 0.5). Mass splitting and the average mass shift tend to zero for y — > 0. 
In this case only the pure gluonic condensates, which enter the even OPE and are 
numerically suppressed compared to other condensates, have a density dependence. 
Note that these evaluations are, at best, for a rough orientation, as mass terms 
oc m s have been neglected. The too low vacuum mass of 1.91 GeV compared to 
the experimental value rri£, s = 1.968 GeV is an indication for some importance of 
strange quark mass terms. Such mass terms oc m s have been accounted for in [21] for 
the vacuum case. The complete in-medium OPE and sum rule evaluation deserves 
separate investigations, as m s introduces a second mass scale. 



V. EVALUATION FOR B AND B MESONS 

We turn now to B and B mesons. The corresponding current operators are 
j B + = or ] B o = ibj^d. The antiparticles correspond to ] B - = ] B+ = iuj^b or 

j-gO = j^o = idr/zb. The above equations and, in particular, the OPE are applied 
with the replacements m c — > m& and m B ± — > m T in order to take into account 
the distinct heavy-light structure compared to the D meson case. The Borel curves 
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0.15 



FIG. 3: Am (left) and m (right) for D s and D s at Sq = 8.0 GeV 2 and for y = 0.5 (solid), 
y = 0.25 (dotted), y = (dashed). 

m±(M) display, analog to the case of open charm, pronounced minima at a Borel 
mass of about 1.7 GeV. We utilize again the first evaluation strategy. Numerical 
results are exhibited in Fig. |4| We employ Sq = 40 GeV 2 and obtain m « 5.33 
GeV for the vacuum mass. One observes a mass splitting of 2 Am ^ — 130 MeV at 
n = 0.15 fm~ 3 . The centroid is shifted upwards by about 60 MeV. The splitting of 
the threshold parameters becomes As 2 , ~ —3.4 GeV 2 and the minima of the Borel 
curves m±(M) are shifted from M « 1.67 GeV in vacuum to M « 1.71 GeV at 
n = 0.15 fnT 3 . In case of B, B mesons, the combination m^idd) is expected to have 
numerically an even stronger impact than the term m c (dd) in the charm sector. 
Indeed, the influence of the chiral condensate becomes even larger than that of the 
odd mixed quark-gluon condensate (q^ga&q) at higher densities. The overall pattern 
resembles the results exhibited in Fig. [2j but with shifted mass scale for m. The other 
evaluation strategies yield the same results. Setting (q^gcr&q) = —0.33 GeV 2 n, and, 
hence, (q^D^q) = —0.0585 GeV 2 n, a mass splitting of 2Am w —220 MeV and an 
average mass shift m 45 MeV would be obtained. 

VI. SUMMARY 



In summary we have evaluated the Borel transformed QCD sum rules for pseudo- 
scalar mesons composed of a combination of a light and a heavy quark. The heavy 
quark mass introduces a new scale compared to QCD sum rules in the light quark 
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FIG. 4: Am (left) and m (right) for B and B at Sq = 40 GeV 2 and m b = 4.7 GeV. For 
line codes see Fig. [7J For density dependent thresholds, Sq = Sq(0) ±7n/no GeV 2 has been 
used. 

sector. The evaluation of the sum rules, complete up to mass dimension 5, has 
been performed for D, D and B, B mesons with a glimpse on D s , D s as well. Our 
analysis relies on the often employed pole + continuum ansatz for the hadronic 
spectral function. This is a severe restriction of the generality of the practical use 
of sum rules. In this respect, the extracted parameters refer to this special ansatz 
and should be considered as indicators for changes of the true spectral functions of 
hadrons embedded in cold nuclear matter. Particles and antiparticles are coupled - 
a problem which is faced also for hadrons with conserved quantum numbers in the 
light quark sector [6j EHJ [30]. 

We presented a transparent approximation to highlight the role of the even and 
odd parts of the OPE. Numerically, we find fairly robust mass splittings (for the 
employed set of condensate values) in the iso-doublets, while an assignment of a 
possible mass shift of the centroids is not yet on firm ground. The impact of var- 
ious condensates is discussed, and (q^q), (q^ga&q) and (qq) are identified to drive 
essentially the mass splitting. While (qq) is amplified by the heavy quark mass, it 
enters nevertheless the sum rules beyond the linear density dependence. A concern 
is the sign of (q^ga&q) , vanishing in vacuum, which determines the size of the D — D 
mass splitting. These findings, in particular for D,D, D S ,D S , are of relevance for 
the planned experiments at FAIR. 
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